arXiv:1505.08124vl [cond-mat.mes-hall] 29 May 2015 


Ag and An Atoms Intercalated in Bilayer Heterostructures of 
Transition Metal Dichalcogenides and Graphene 

F. lyikanat/’Q H. Sahin,^’^ R. T. Senger/ and F. M. Peeters^ 

^Department of Physics, Izmir Institute of Technology, 35430, Izmir, Turkey 
^Department of Physics, University of Antwerp, 2610, Antwerp, Belgium 

(Dated: June 1, 2015) 

Abstract 

The diffusive motion of metal nanoparticles Au and Ag on monolayer and between bilayer het¬ 
erostructures of transition metal dichalcogenides and graphene are investigated in the framework 
of density functional theory. We found that the minimum energy barriers for diffusion and the 
possibility of cluster formation depend strongly on both the type of nanoparticle and the type of 
monolayers and bilayers. Moreover, the tendency to form clusters of Ag and Au can be tuned by 
creating various bilayers. Tunability of the diffusion characteristics of adatoms in van der Waals 
heterostructures holds promise for controllable growth of nanostructures. 
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I. INTRODUCTION 


In recent years graphene^ has become one of the most attractive materials due to 
its unique properties such as high-mobility electron transport,-*^ the presence of room- 
temperature quantum Hall effect,- the strong lattice structure^ and the extremely high 
in-plane thermal conductivity.- However, its highly active surface and the lack of a band 
gap in the electronic structure are emerging drawbacks for graphene. Recently, interests 
have now also focused on other two-dimensional systems having honeycomb structures, such 
as graphane,-'^ halogenated graphenes,-”— silicene,— HI-V binary compounds,— and tran¬ 
sition metal dichalcogenides (TMDs).— ~— Recent studies have revealed that among various 
monolayer structures especially TMDs are quite promising materials for electronics and 
optoelectronics applications. 

Bulk TMDs have a number of exceptional properties such as superconductivity of TaS2 
and NbSe2, Mott transition in lT-TaS2 and the presence of charge density wave in TiSe2.— 

It was also reported that the electrical and optical properties of TMDs are dramatically 
altered with the number of layers.—”— Although bulk hexagonal TMDs possess an indi¬ 
rect band gap, mono-layer TMDs exhibit a direct band gap which is crucial for optoelec¬ 
tronic devices, sensors and catalysts. In addition, n-type and p-type held-effect-transistors 
(FETs) based on monolayer and multilayer TMDs have been investigated.—”— It was also 
reported that many monolayer 2 D crystals are reactive and segregation may occur easily.— 
Therefore, the investigation of bilayer structures which are chemically more stable than 
monolayer structures is of vital signihcance. Recent studies have shown that synthesis 
of heterostructures made of combinations of different TMD single layers, graphene, fluo- 
rographene and hexagonal-BN (hBN) is experimentally achievable.—”— Since TMDs and 
other two-dimensional structures have a lot of diverse monolayer structures, when they are 
combined together they are expected to exhibit very different properties.—*^ 

Since the intercalation and migration of foreign atoms is inevitable during the formation 
of such lamellar materials and heterostructures, the investigation of the diffusion characteris¬ 
tics of various impurities is essential. Early studies revealed that alkali-metal doping of bulk 
TiSe2 and M0S2 can be utilized as an efficient way to tune the Fermi level.— The electrical 
conductivity of M0S2 can be altered by substitutional doping.— > 2 ^ Furthermore, decoration 
of the surfaces of few layer TMDs by metal nanoparticles like An, Ag and Ft may provide p- 
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and n-type doping.—^— The metal-atom adsorbed TMDs find their application in various ar¬ 
eas including energy storage,— photonics,—*^ biosensing,— and catalysis.— In a recent study, 
it has been demonstrated that, the presence of various impurities at the interface between 
MoS 2 /graphene/hBN and WS 2 /graphene/hBN heterostructures may modify the mobility of 
charge carriers.— It was also reported that contamination and migration of various molecules 
are inevitable during the formation of graphene based heterostructures and trapped hydro¬ 
carbons segregate into isolated pockets, leaving the rest of the interface atomically clean.— 
In addition, it was found that, attached metal nanoparticles on TMDs/Graphene stacks can 
be suitable for enhanced optoelectronic properties.—*^ 

Despite some recent studies on adatom adsorption on various TMDs, intercalation and 
migration of foreign atoms in heterostructures have not been investigated. In this study, 
using density functional theory based electronic structure method, we investigate the diffu¬ 
sion characteristics of heavy metal atoms (An and Ag) on monolayers and intercalated in 
such bilayer heterostructures. The paper is organized as follows: In Sec. [TTl we give details 
of our computational methodology. In Sec. Illll the energetics of the metal atoms Ag and An 
on monolayers of graphene and TMDs are presented. In Sec. IIVI diffusion characteristics of 
those metal atoms inside bilayer heterostructures are shown and in Sec. |V] we summarize 
and conclude our results. 

II. COMPUTATIONAL METHODOLOGY 

To determine ground state atomic structures and migration characteristics of monolayers 
and their bilayer heterostructures, first-principles calculations were performed using den¬ 
sity functional theory (DFT) with a plane-wave basis set as implemented in the Vienna ab 
initio simulation package (VASP).— For the exchange correlation function generalized gra¬ 
dient approximation (GGA) of Perdew, Burke and Ernzerhof^ was used together with the 
van der Waals correction.— Spin-unpolarized calculations were carried out using projector- 
augmented-wave potentials (PAW). The plane-wave basis set with kinetic energy cutoff of 
500 eV was used. 

A 4x4 and 3x3 hexagonal supercells of single layer and bilayer structures are employed 
to model diffusion paths of metal atoms, respectively. The k-point samples were 3x3x1 for 
these supercells. It is calculated that a 12 A of vacuum space is enough to hinder long-range 
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dispersion forces between two adjacent images in the supercell. Lattice constants and total 
energies were computed with the conjugate gradient method, where atomic forces and total 
energies were minimized. The convergence criterion of our calculations for ionic relaxations 
is 10“^ eV between two consecutive steps. The maximum Hellmann-Feynman forces acting 
on each atom were reduced to a value of less than 10“^ eV/A. All components of the external 
pressure in the unit cell was held below 1 kbar. For the electronic density of states Gaussian 
smearing was used with a broadening of 0.1 eV. The non-local correlation energies were 
determined by employing density functional theory plus the long-range dispersion correction 
(DFT+D2) method.-^ The values of 24.670, 5.570, 12.640, 1.750, 24.670, 40.620, and 40.620 
are used as the Cg coefficient of Mo, S, Se, C, Ag, Au, and W atoms, respectively. On the 
other hand, values of 1.639, 1.683, 1.771, 1.452, 1.639, 1.772, and 1.772 are used for the vdW 
radius of Mo, S, Se, C, Ag, Au, and W atoms, respectively. These Cg coefficients and vdW 
radii were determined from previous studies.— It is essential that long-range dispersion 
correction to the interlayer force is taken into account in order to obtain reliable layer-layer 
distances and electronic properties of the heterobilayers. 

To study the adsorption and diffusion of metal atoms in these systems, total energies 
of monolayer TMDs (graphene) with adatom were calculated for 19 (12) different points 
which include the high symmetry points (H, M, B, and T). Total energies of bilayers with 
adatom were calculated for only 4 high symmetry points. In these different points, the 
adatom-surface distance was fully relaxed while the position of the adatoms parallel to 
the surface was kept hxed. To obtain accurate diffusion characteristics of heavy atoms 
on monolayer TMDs and graphene, calculations were performed using 4x4 supercells. In 
these calculations, first and second nearest neighbors of the adatom were fully relaxed while 
all the rest was kept fixed. On the other hand, in the bilayer calculations, all atoms of 
bilayers were free to move in all directions. Binding energies were calculated for the most 
favorable adsorption sites. These binding energies were calculated from the expression Eb = 
EMonolayer + Ea — EMonolayer+A, where Efi is the binding energy of the metal atoms on the 
TMDs or graphene, ^Monolayer is the energy of monolayer TMDs or graphene, is the 
energy of metal atoms, EMonoiayer+A is the total energy of the metal atom-monolayer system. 

We have chosen the following convention to define the diffusion energy barrier. The 
binding energy of an adatom was calculated at all high-symmetry points of the TMDs and 
the graphene surface. Since adatoms follow the lowest energy path, the difference of the 
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energies between the most favorable site and the second most favorable site was considered 
as the diffusion energy barrier. 

In order to obtain the correct value of the charge transferred between the metal atoms, 
the TMDs and the graphene, Bader charge analysis was performed. Rather than electronic 
orbitals in Bader methodology, charge partitioning is based on electronic charge density and 
therefore it is highly efficient, it scales linearly with the number of grid points, and is more 
robust than other partitioning schemes.—^— 


III. DIFFUSION OF HEAVY ATOMS ON MONOLAYERS OF TMDS AND 
GRAPHENE 

So far, experiments have revealed that similar to graphene many of the TMDs (M 0 S 2 , 
MoSe 2 , WS 2 , WSe 2 , ReS 2 ) have a lamellar crystal structure where the layers are held together 
by weak van der Waals forces, while intra-layer metal-chalcogen bonds have strong covalent 
character. Due to the weakly (vdW) bonded layered structure of graphene and TMDs 
either diffusion of foreign atoms on their surface or easy intercalation into the layers may 
take place. Moreover, these foreign atoms may provide new functionalities to the TMDs and 
to the heterostructures made of various layered materials. Therefore, understanding how 
foreign atoms adsorb, migrate and intercalate are of highly importance. 

First we calculate the binding and migration of Ag and An atoms on the surface of 
monolayer TMDs and graphene. For a hexagonal primitive unitcell of M 0 S 2 , MoSe 2 , WS 2 
and graphene the lattice parameters are calculated to be 3.190, 3.316, 3.179 and 2.468 A, 
respectively. For calculations of adsorption and migration of metal atoms 4x4 supercells 
were used in order to limit Au-Au (and Ag-Ag) interaction between adjacent cells. As 
shown in Fig. [T]for a metal atom adsorbed on TMDs and graphene, there are four possible 
adsorption sites: top of carbon or chalcogen atom (T), bridge site on C-C or M-X bond (B), 
on top of the center of hollow site (H) and on top of metal atom Mo or W (M). 

Diffusion and mobility of impurities on a two-dimensional crystal surface can be described 
by quantities such as activation energy {Ea) and jump probability (P) from the binding site 
over the lowest-barrier path at room temperature. The P value is a measure of the possibility 
of propagation by overcoming the energy barrier among the possible adsorption sites. Jump 
probability from one lattice site to another one can be calculated by using the formula 
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P ^ g-Ea/fcsT ks is the Boltzmann constant and it increases exponentially with 

increasing temperature. Here Ea is the activation energy which is equal to the difference of 
the energy of the two lowest energy states. 

Our calculated adsorption sites, binding energies, vertical adsorption position, charge 
transfer, energy barrier and jump probability (P) values are listed in Table [H We see that 
the bonding of an Ag atom to graphene occurs at the B site and the H site is the least 
favorable adsorption site. In accordance with the previous DFT study, we found that the 
binding energy of Au on graphene is almost twice the binding energy of Ag on graphene. The 
Ag atoms lose (while the Au atoms gain) a small amount of charge when they bind to the 
graphene surface.— On the surface of TMDs, the most favorable bonding site for Ag atoms 
is the H site and the next largest adsorption energy is found for the M site. It is reasonable 
to assume that metal atoms diffuse through these two favorable adsorption sites. Therefore 
the energy difference between these two lowest-energy sites can be regarded as the diffusion 
barrier. Since the energy barrier for Ag between H and M sites is high (~50 meV) diffusion 
through these symmetry points may not occur at low temperatures. It appears from Fig. [2] 
that Ag atoms on graphene migrate through the B and T sites with almost zero (~2 meV) 
energy barrier, while migration on TMDs (M 0 S 2 , MoSe 2 and WS 2 ) occurs through H and 
M sites. As shown in Table [H room temperature P values of Ag atoms are 0.129, 0.182, 
0.072 and 0.926 for M 0 S 2 , MoSe 2 , WS 2 and graphene, respectively. Thus diffusion of the Ag 
atoms on graphene occurs much more easily than on TMDs. 

Binding energy of a single Au atom on graphene and TMDs is found to be ~300 meV 
larger than that of Ag. Interestingly, on MoSe 2 , Au atoms prefer bonding on the H site like 
Ag atoms, while the top of a sulfur atom (T site) is the most preferable site on M 0 S 2 and 
WS 2 layers. The distinctive behavior of the Au atoms can be explained by the distinctive 
characteristic of the S-Au bond. Although the chemical properties of chalcogen atoms S 
and Se are quite similar, as shown by Yee et ah,— the strength of the S-Au bond is slightly 
weaker than the Se-Au bond. Our calculations reveal that the binding energy of an Au 
atom on MoSe 2 surface is higher than the binding energy of an Au atom on M 0 S 2 and WS 2 
surfaces. Interestingly, strong Au-S bonds require bonding on T site while the Au-Se bond 
results in H-site adsorption. As shown in FigjSl as in the Ag adsorption case, M 0 S 2 and WS 2 
surfaces have almost identical characteristics for Au adsorption at all high symmetry points. 
In addition, the migration barrier seen by Au atoms is ~30 meV larger than the energy 
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TABLE I: Several adsorption characteristics of adatoms on monolayer substrates. Binding site 
and A/9 denote the energetically most favorable position and change in the charge of adatom, 
respectively. Binding energy, height (relative to surface), energy barrier (at high symmetry sites, 
relative to the binding site) and jump probability (P) of the most favorable sites of Au and Ag 
atoms on the TMDs and graphene surface are listed. 



Ag/MoS2 

Ag/MoSe 2 

Ag/WS2 

Ag/Graphene 

Binding Site 

H 

H 

H 

B 

Binding energy (eV) 

0.943 

0.928 

0.845 

0.360 

Height (A) 

1.97 

2.07 

2.04 

2.57 

^P = Pf - P* (e) 

-0.30 

-0.18 

-0.23 

-0.16 

Energy Barrier (meV) 53(M), 146(B), 228(T) 44(M), 139(B), 241(T) 68(M), 111(B), 181(T) 

2(T), 57(H) 

Jump probability (P) 

0.129 

0.182 

0.072 

0.926 


AU/M 0 S 2 

Au/MoSe 2 

Au/WS2 

An/Graphene 

Binding Site 

T 

H 

T 

T 

Binding energy (eV) 

1.161 

1.232 

1.146 

0.633 

Height (A) 

2.26 

2.03 

2.27 

2.56 

Ap = p/ - Pi (e) 

0.01 

0.08 

0.06 

0.22 

Energy Barrier (meV) 

136(M), 23(H), 46(B) 

103(M), 72(B), 68(T) 

180(M), 43(H), 55(B) 

65(H), 26(B) 

Jump probability (P) 

0.411 

0.072 

0.189 

0.366 


barrier of Ag atoms on graphene. Therefore, compared to Ag on graphene, the smaller jump 
probability of Au atoms on graphene implies relatively slow and smaller nucleation of Au 
clusters on graphene which agrees with experimental results.— 

Calculated room temperature P values of an Au atom on M 0 S 2 , MoSe 2 , WS 2 and 
graphene surfaces are 0.411, 0.072, 0.189 and 0.366, respectively. Au atom has the maxi¬ 
mum P value on a M 0 S 2 surface, whereas it has the minimum P value on a MoSe 2 surface. 
P values of Au atoms on a M 0 S 2 surface is higher than that of Ag atoms. This result is 
consistent with the previous experimental study where they found that the cluster diameter 
of Au and Ag atoms on M 0 S 2 were found to be 14 and 6 nm, respectively.— Other exper¬ 
imental work showed that cluster diameters of Au and Ag atoms on WS 2 are 19-20 and 5 
nm, respectively.— In agreement with this experimental study, our calculations show that 
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the P value of the Au atoms on a WS 2 surface is higher than that of Ag atoms. 


IV. DIFFUSION OF HEAVY ATOMS INTERCALATED BETWEEN VAN DER 
WAALS BILAYERS 

The vertical stacking of graphene and other two-dimensional atomic crystals allows for the 
combination of different electronic properties. Recent studies have shown that, electronic and 
optical properties of TMDs can be altered dramatically by forming such heterostructures.— 

It appears that the intercalation and contamination of foreign atoms and functional groups 
at the interface of these heterostructures are unavoidable. As shown by recent experimental 
study, heavy atoms such as Au and Ag are quite mobile on TMD surfaces and they form 
clusters.— Therefore, understanding the diffusion characteristics of heavy atoms at the in¬ 
terface of these heterostructures is of importance for the ongoing research on heterostructure 
devices. 

In this section, we investigate the diffusion characteristics of Au and Ag atoms intercalated 
in MoS 2 /Graphene, M 0 S 2 /M 0 S 2 , MoS 2 /MoSe 2 and M 0 S 2 /WS 2 heterostructures. 

A. Diffusion through MoS2/Graphene Bilayer Heterostructure 

First, we start with the heterostructure MoS 2 /Graphene. However, the determination of 
the most favorable atomic structure of MoS 2 /Graphene heterostructure is more complicated 
due to the lattice mismatch. Our calculations showed that bare DFT calculations are not 
capable of hnding the ground state ordering of graphene placed on M 0 S 2 due to the presence 
of many local minima corresponding to metastable states. It was seen that low-temperature 
(100 K) molecular dynamic calculations are able to avoid local minima allowing to obtain 
the relaxed geometric structure of MoS 2 /Graphene as shown in Fig. [I](c). Due to the lack of 
symmetry in the MoS 2 /Graphene structure, twenty-seven inequivalent adsorption sites are 
considered for each Ag (and Au) atom adsorption. Table [ITl lists the intercalation energies 
obtained by using a 3x3 unit cell of the most favorable site for each configuration. Here the 
intercalation energy is dehned by the expression 

Pint pHetero+Ag(Au) pRetero pAg{Au)i 


( 1 ) 


where Eint is the intercalation energy of Ag(Au) atom, EHetero+Ag{Au) is the energy of the 
metal-atom heterobilayer system, Euetero is the energy of heterobilayers and EaqIAu) is the 
energy of the adatom. Due to the large binding energies of adatoms on M 0 S 2 , all the 
equilibrium geometries are ruled by M 0 S 2 except for the MoS 2 /Au/MoSe 2 structure where 
Au’s binding is larger with the MoSe 2 layer. 

Our calculations revealed that binding energies of Ag and Au atoms on TMD substrates 
are larger than that on graphene and Au atoms are more stable than Ag atoms on these 
monolayers. For the simulation of the MoS 2 /Graphene heterostructure we considered 4x4 
and 3x3 unit cells for graphene and M 0 S 2 layers, respectively. After full optimization, the 
lattice parameter of the heterostructure is calculated to be 9.791 A. For this supercell the 
maximum lattice mismatch is ~3%. 

The adsorption energies and the equilibrium geometries were calculated by considering 
27 different adsorption sites for each Ag and Au atoms at the MoS 2 /Graphene interface. 
It is seen that the most favorable site of Ag atoms on M 0 S 2 is changed from H to M site 
upon the addition of the graphene layer. Energy barrier of Ag atoms is enhanced from 2 to 
43 meV with graphene layer. Although the energy barriers seen by Au atom on M 0 S 2 are 
significantly enhanced by the presence of the graphene layer, the T site of M 0 S 2 remains as 
the most favorable adsorption site of Au. 

Furthermore, we show 2D plots of the energy barriers seen by Ag and Au atoms at 
the MoS 2 /Graphene interface in Fig. IH It is clear that Ag and Au atoms have entirely 
different diffusion characteristics at the M 0 S 2 /Graphene interface: (i) diffusion of Ag atoms 
are much easier while Au atoms may diffuse at high temperatures, (ii) diffusion of an Ag 
atom occurs through the bridge sites of the top-lying graphene layer, (hi) diffusion of Au 
atom may occur from one top-sulfur site to another one through the top site of Mo atoms. 
Here, similar characteristic diffusion behavior of heavy atoms can be expected for bilayer 
heterostructures of WS 2 /Graphene as well. 

B. Diffusion through MX2/MX2 Bilayer Heterostructures 

Next we investigate the diffusion and energy barriers which are seen by Ag and Au atoms 
sandwiched in between single layers of TMDs. There are two possible stacking types for 
bilayer TMDs; AA and AB stacking. In the AA stacking, the chalcogen atom (S or Se) in 
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TABLE II: Several adsorption characteristic of adatoms between bilayer heterostructures. Binding 
site and Ap denotes the energetically most favorable position and change in the charge of adatom, 
respectively. LL Distance gives the interlayer distance. Intercalation energy, height above M 0 S 2 , 
energy barrier (at high symmetry sites, relative to the binding site) and jump probability (P) of 
the most favorable sites of Au and Ag atoms between bilayer heterostructures are listed. 



MoS2/Ag/MoS2 

M 0 S 2 / Ag/MoSe 2 

MoS2/Ag/WS2 

MoS 2 /Ag/Graphene 

Binding Site 

H 

H 

H 

M 

Intercalation Energy (eV) 

-1.461 

-1.395 

-1.319 

-1.344 

Height above M 0 S 2 (A) 

1.85 

1.71 

1.85 

1.61 

Ap = p/ - p^ (e) 

-0.42 

-0.41 

-0.41 

-0.47 

LL Distance (pristine, with Ag) 

3.08, 3.70 

3.10, 3.60 

3.09, 3.71 

3.32, 3.86 

Energy Barrier (meV) 

222(M), 282(B) 

206(M), 308(B) 

243(M), 295(B) 

43(H), 196(T) 


222(T) 

274(T) 

340(T) 


Jump probability (P) 

2x10-4 

3.5x10-4 

8x10-5 

0.189 


M 0 S 2 /AU/M 0 S 2 

MoS2/Au/MoSe2 

M 0 S 2 /AU/WS 2 

M 0 S 2 / Au/Graphene 

Binding Site 

B 

H 

H 

T 

Intercalation Energy (eV) 

-1.645 

-1.562 

-1.504 

-1.541 

Height above M 0 S 2 (A) 

2.07 

1.79 

1.85 

2.07 

^P = Pf - P* (e) 

-0.15 

-0.13 

-0.18 

-0.21 

LL Distance (pristine, with Au) 

3.08, 4.12 

3.10, 3.72 

3.09, 3.74 

3.32, 4.35 

Energy Barrier (meV) 

12(H), 96(M), 95(T) 118(M), 75(B), 103(T) 78(M), 24(B), 164(T) 

468(H) 326(M) 

Jump probability (P) 

0.628 

0.055 

0.395 

3x10-® 


one layer is on top of the chalcogen atom (S or Se) of the other layer. On the other hand, in 
AB stacking, chalcogen atom (S or Se) in one layer is on top of the metal atoms (Mo or W) 
of the other layer. Our calculations show that for the considered bilayer heterostructures of 
TMDs, the AB stacking is the ground state geometry. 

M 0 S 2 /M 0 S 2 Bilayers: The lattice parameter of AB-stacked bilayer M 0 S 2 /M 0 S 2 structure 
was found to be almost the same as of monolayer M 0 S 2 . It is seen that the Ag-MoS 2 distance 
decreases due to the presence of the upper M 0 S 2 layer. The layer-layer distance of bilayer 
M 0 S 2 increases from 3.08 to 3.70 A upon single Ag intercalation. On the other hand. 
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single An atom causes a larger layer separation (from 3.08 to 4.12 A) than Ag atoms. Our 
calculations revealed that the most preferable site of the Ag atom remains the H-site, even 
after the addition of the second layer. Our charge analysis showed that metal atoms tend 
to donate a larger amount of charge when they are intercalated in bilayers. The amount 
of charge transfer becomes 0.42e and O.lSe for Ag and An atoms, respectively. However, it 
is worthwhile to note that the preferable adsorption site for An atoms changes from T to 
B site. An atoms tend to bind with two S atoms that are in two different layers. For this 
reason, the B-site is the most favorable site. 

The minimum diffusion barrier of Ag atom is increased by ~ 170 meV after the addition 
of the second layer. Interestingly diffusion barriers of the Ag atom between bilayer M 0 S 2 
from H to M and H to T sites are the same. The minimum diffusion barrier of An atom is 
reduced by half when a second layer is added. These results show that the diffusion behavior 
of Ag and An atoms between bilayer M 0 S 2 are completely different. Ag atoms are not likely 
to form clusters inside the bilayer M 0 S 2 . On the other hand. An atoms may form clusters 
between two M 0 S 2 layers and the clusterization weight of them is larger than that on single 
layer of M 0 S 2 . 

MoS 2 /MoSe 2 Bilayers: For the geometry optimization of the MoS 2 /MoSe 2 bilayer struc¬ 
ture we performed full atomic and lattice parameter relaxation and therefore internal strain 
due to lattice mismatch is lowered to less than 1 kB. The lattice parameter of MoS 2 /MoSe 2 is 
calculated to be 3.252 A. For this lattice parameter, the maximum lattice mismatch is ~2% 
which is in the range of experimentally available strain values. When An and Ag atoms are 
intercalated between MoS 2 /MoSe 2 , they show the same trend; their height above the M 0 S 2 
surface are fairly close (1.7 A for Ag and 1.8 A for An) and the layer-layer distance increases 
from 3.1 to 3.6 and 3.7 A for Ag and An, respectively. H-site is the most favorable site for 
Ag and An atoms when they are located between MoS 2 /MoSe 2 bilayer. When Ag atom is 
inserted in MoS 2 /MoSe 2 , it donates 0.41e to neighboring TMD layers. It appears that the 
majority of the charges are transfered to the M 0 S 2 layer. Similarly the An atom loses 0.13e 
charge when inserted between MoS 2 /MoSe 2 . The minimum diffusion barrier of Ag atom is 
increased by ~ 150 meV when the second layer is present. But the second favorable site 
does not change with the addition of the MoSe 2 layer (M). The minimum diffusion barrier 
of An atom increases from 23 to 75 meV with MoSe 2 addition. Second favorable site of An 
atom changes from H to B site when MoSe 2 layer is added. Similar to M 0 S 2 /M 0 S 2 , the jump 
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probability of Ag atom becomes zero with the addition of MoSe 2 as a second layer. The P 
value of An atom decreases from 0.411 to 0.055 with the addition of MoSe 2 layer. Therefore, 
the formation of Ag cluster between MoS 2 /MoSe 2 bilayer heterostructure is hindered by 
the presence of the second layer. The nucleation of An clusters may happen with a slower 
nucleation rate in MoS 2 /MoSe 2 as compared to monolayer M 0 S 2 . 

M 0 S 2 /WS 2 Bilayers: The lattice mismatch between M 0 S 2 and WS 2 is negligible and the 
lattice parameter of the bilayer structure was calculated to be 3.182 A. Both An and Ag 
atoms prefer to be bonded on the H site 1.85 A away from the surface of M 0 S 2 . When 
Ag and An atoms are inserted between M 0 S 2 /WS 2 , the interlayer distances increase from 
3.09 to 3.71 and 3.74 A, respectively. We found that the Ag atom donates 0.41 e to the 
surrounding M 0 S 2 /WS 2 bilayer structure whereas the amount of charge transfer is 0.18 e for 
An. The minimum diffusion barrier of the Ag atom increases ~ 190 meV with the addition 
of WS 2 layer, but it does not change the second favorable site of Ag atom (M). It is seen that 
the addition of the WS 2 layer does not have any signihcant effect on the minimum diffusion 
barrier of the An atom. 

In this case, the change of the P value of the Ag atom shows the same trend with bilayer 
M 0 S 2 and MoS 2 /MoSe 2 . In these cases, the P value of Ag atom becomes zero with the 
addition of the second layer (WS 2 in this case). On the other hand, addition of a WS 2 
layer does not change the P value of the An atom between M 0 S 2 /WS 2 , significantly. As a 
result, Ag atoms are not likely to form clusters when intercalated in M 0 S 2 /WS 2 , while the 
clusterization probability of An atoms does almost not change with WS 2 layer addition. 

V. CONCLUSION 

We found that the binding energy of Ag and An atoms is maximum when they are placed 
on M 0 S 2 and MoSe 2 , respectively and the binding energy of these atoms on graphene is lower 
than that on TMDs. Ag atoms prefer to bind at the H site on TMDs, while the H and the 
T sites are most favorable adsorption sites for An atoms. Just as for the Ag atom. An atom 
binds to H site on MoSe 2 . However, An atom binds to T site on M 0 S 2 and WS 2 . When 
Ag atom adsorbs on these monolayers, it tends to donate charge, whereas An atom tends 
to receive charge from these monolayers. Calculated P values revealed that, clusterization 
tendency of Ag and An are the largest on graphene and M 0 S 2 layers, respectively. 
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Addition of a second layer has a considerable effect on the diffusion phenomena of Ag 
and An atoms. Ag atom binds to H site of all TMDs heterostructures, whereas it binds 
to M site of the MoS 2 /Graphene heterostructure. When an An atom is inserted between 
MoS 2 /MoSe 2 and M 0 S 2 /WS 2 heterostructures, H site is the energetically most favorable 
site. B and T sites are the most favorable sites for the An atom when it is placed between 
M 0 S 2 /M 0 S 2 and MoS 2 /Graphene heterostructures. The clusterization probability of Ag 
atom becomes zero when TMDs are added as a second layer. On the other hand, the P 
value of Ag atom increases with the addition of a graphene layer. Addition of a graphene 
layer hinders diffusion of An atom. WS 2 addition as a second layer, does not make a 
signihcant effect on the cluster formation probability of An atom. Addition of MoSe 2 as a 
second layer reduces the clusterization of An, whereas the addition of a second M 0 S 2 layer 
increases it. Our calculated intercalation energies showed that, when An and Ag atoms are 
inserted between M 0 S 2 /Graphene, M 0 S 2 /M 0 S 2 , MoS 2 /MoSe 2 or M 0 S 2 /WS 2 , they are ener¬ 
getically more stable as compared to the case where Ag and An atoms are adsorbed on the 
monolayer substrates. The energy barrier for An atoms is lower than the one for Ag atoms, 
when they are inserted between TMD heterostructures. On the other hand, the energy 
barrier of Ag atoms is lower than the energy barrier of An atoms when they are inserted 
between MoS 2 /Graphene. Tunability of diffusion barriers seen by adatoms in van der Waals 
heterostructures makes these materials potential hosts for controlled growth of nanoclusters. 
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Figl: (Color online) Possible adsorption sites on (a) TMD and (b) Graphene. Top (left) 
and side (right) views of bilayer heterostructures of (c) TMD/Graphene and (d) TMD/TMD. 

Fig2: (Color online) Diffusion paths of Ag on Graphene, M 0 S 2 , MoSe 2 , and WS 2 . (see 
Fig {T] for sites) 

Fig3: (Color online) Diffusion paths of An on Graphene, M 0 S 2 , MoSe 2 , and WS 2 . (see 
FiglUfor sites) 

Fig4: (Golor online) Gontour plots of the energy barriers (in meV units) seen by Ag 
(upper panel) and An (lower panel) atoms intercalated in M 0 S 2 /Graphene heterostructure. 
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